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The  investigation  of  drag-characteristic  behavior  of  the  wing  in 
steady  flow  is  omitted  entirely  from  this  paper^  because  it  is  not  a 
difficult  matter  for  a  designer  to  calculate  the  drag  characteristics 
with  the  information  contained  in  this  paper  and  available  solutions 
of  the  steady  case. 

The  present  investigation  was  conducted  at  the  Johns  Hopkins 
University  under  the  sponsorship  and  with  the  financial  assistance  of 
the  National  Advisory  Committee  for  Aeronautics.  The  author  would  like 
to  express  his  appreciation  to  Misses  Vivian  O'Brien  and  Patricia 
darken  for  their  assistance  in  carrying  out  the  project. 


TRANSIENT  REACTION  OF  AN  AIRFOIL  DUE  TO  UNIT -STEP  CHANGE 
OF  FLAPPING  AND  VERTICAL  GUST 


As  shown  in  references  1  and  6,  the  linearized  partial  differential 
equation  of  the  Irrotational  flow  of  a  compressible  nonviscous  fluid  is 


where  0  is  the  disturbance  potential  of  the  flow,  the  x-axis  is 
positive  in  the  opposite  direction  of  the  flight,  and  the  y-axis  is 
positive  in  the  upward  direction.  The  origin  is  attached  to  the  air- 
edge.  The  details  of  the  notations  are  given  in 
appendix  A. ^  In  the  supersonic  case,  the  disturbance  potential  of  such 
an  airfoil  is  equivalent  to  that  due  to  a  source  sheet  at  the  y  =  0 
plane  with  time -dependent  strength  and  is 


0(x,y,t)  =  -  — 

jt 


1 

(m2  -  1)1/2 


v( g^+0,t  -  t) 

]j(T  -  Ti)  (t2  -  t) 


dr  d| 


(2) 


where  v( |,+0,t  -  t)  Is  the  vertical  velocity  component  on  the  wing 
surface  at  the  point  (l,+0)  at  an  earlier  time  t  -  t.  The  term 
v( l,+0,t  -  t) 

—  ^  ,  is  just  the  so\irce-sheet  strength  per  unit  area.  In 

(m2  -  1)  / 

addition, 


k 
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where 


If  the  point  of  interest  is  on  the  lower  side  of  the  wing  surface 
(x,-0^t)  equation  (2)  reduces  to 


(5) 


A 


With  the  above  eq\iation  the  transient  case  due  to  sudden  change  of 
angle  of  attack  of  the  airfoil  about  any  arbitrary  point  has  been  treated 
in  reference  1.  The  important  results  of  that  paper  are  given  in  table  1 
for  use  in  the  present  development.  Two  additional  interesting  cases 
will  be  treated  as  follows. 


Vertical  Flapping 

At  time  tQ,  the  airfoil  suddenly  descends  with  constant  velocity 
The  unit-step  downwash  at  any  point  on  the  wing  can  be  expressed  in  terms 
of  Dirichlet’s  integral  as 


v(x^-O^t)  =  fiiCt)  =  fi  — 


sin  03 


(t  -  to) 


dm  +  i 
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where 


u(t  -  to)  "  (f  -  to') 


-  1 


(m^  -  ^ 


to’  = 


and  jQ(a)\)  is  the  Bessel  function  of  the  first  kind  of  the  zero  order. 

^h  ^nh  ^Lh 

Table  2  gives  the  detailed  characteristics  of 

aji  ah  ah 

Cj4h  , 

and  n —  =  as  functions  of  time.  Figures  2  to  4  shov  these  charac- 

% 

teristics.  The  detailed  methods  of  calculation  are  quite  similar  to 
those  given  in  reference  1,  and  are  omitted  in  this  paper. 


In  figures  2  to  4^  a  few  important  features  may  be  explained.  Within 


the  zone  III  { 0  ^ 


M  \  ^ph 

- L?  f-h  =  is  inversely  proportional 

M  +  1/  % 


to  Mach  number  and  independent  of  time  and  the  location  of  x* .  But  in 


time  zone  I 


1 1  t  *  1  C  h 

is  inversely  proportional  to 

M  -  1  x'  J  % 


^2  ^1  which  corresponds  to  the  steady  case  of  Ackeret’ s  result.  In 


zone  II 


M  <  ~  '*^o'  <  _ M 


M  +  1 


— — — I,  varies  not  with  x'  or  t'  -  to' 

M  -  1/^  ail  ^ 


alone,  but  with  the  conical  coordinate 


It  is  a  complicated 


function  of  M  and 


Thus  the  present  solution  is  similar  to 


Busemann’ s  conical  flow.  So  is  the  behavior  of  Fh  z:: —  against  t’  -  tQ* 

ah 
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As  far  as  ^  is  concerned,  it  varies  both  with  1/m  and  t'  -  t^' 

h  0 

in  zone  I.  In  general,  fj^,  Fj^^,  and  Gj,  increase  with  decreasing  Mach 
"t  ^  "t  ^ 

number,  if  or  t'  -  to'  is  kept  constant. 


Vertical  Gust 

OTic:i-  nr  airfoil  begins  to  encounter  a  uniform  vertical 

gust  of  velocity  g,  such  unit-step  downwash  g  at  any  point  on  the 
airfoil  can  also  be  expressed  in  terms  of  Dirichlet' s  integral 


v(x,-0,t)  =  g(t)  =  g  i 


Sin  CJD 


(t  -  to  - 


dco  +  i 
2 


Figure  1(b)  shows 


a-(t)  as  a  function  of  both  t  and  x 


r.n+  equation  (9)  into  equation  (5),  the  .disturbance 

potential  0g  due  to  vertical  gust  can  be  written  as 


(10) 
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or,  after  integration,  0g  has  to  be  also  expressed  in  three  zones. 


It  is  interesting  to  note  that  0h  0g  similar  in  fom.  The 

Fourier  integral  method  of  Bessel  functions  used  will  take  care  of  the 
three  zones  automatically  without  considering  each  zone  individually. 
This  is  one  of  the  main  advantages  of  the  present  method. 


Table  3  shows  the  detailed  characteristics  of 


^  ^  F  and  ^  =  Gg  aS  functions  of  time.  Figures  5  to  7  show 

fg;  and  Gg  as  functions  of  time.  The  detailed  calculations  are 

also  omitted. 


Pitching 

For  convenience^  the  essential  results  of  reference  1  for  the 
pitching  case  are  given  in  table  1.  Figures  8  to  10  show  some  charac¬ 
teristics  of  this  case. 

In  the  case  where  xq*  =  0  or  the  airfoil  rotates  about  the  leading 


'dge  at  a  constant  rate  ^a-|_  -  clq^  starting  at  t^  - 


X* (di  -  6^) 


varies  linearly  with  ^  but  inversely  with  M  in  zone  III 

X* 


t’  -  tQ’  ^ 

x'  ""  M  +  1 


if  the  angle  of  attack  is  small  at  any  instant 
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In  zone  I 


/  M  ^  -  ^0'\ 

\M  -  1  ■  x'  J 


it  also  varies  linearly  with 


f  -  to' 


^  but 


inversely  with  -1.  If  M>1,5^ 


'pa 


x*(ai  -  oo) 


also  varies  approxi- 


t'  -  t 


mately  linearly  with 


0 


and  is  continuous,  in  value  with  zones  I 


X 


and  III  although  the  exact  expression  as  shown  in  table  1  is  quite  com¬ 
plicated.  Of  course,  as  the  Mach  number  is  near  to  1,  this  linear 
approximation  becomes  worse.  It  seems  better  to  use  the  exact  solution 

in  the  table.  Similar  statements  can  be  madt  about  F-  =  — 


(ai  -  Oo) 


and  Gi  = 


'Md 


(di  -  do) 

If  Xo’  i  0  some  additional  contributions  have  to  be  considered  as 


fd  = 


'pa 


X‘  -  do) 


=  f .  »  -  Af*  * 
■^a  ^a 


(12) 


F-  = 


'La 


a  /-r  TT  \  -"a 


=  Fa’  •  AFi’ 


(13) 


%d 


Gx  =  - - =  G^,'  - 


"  (di-do) 


(1^^) 


^d'  ^d’ 

where  — —  =  f>,,  — — 


^d* 


X/ 


X/ 


=  Fj^,  and  -  =  are  given  in  table  2  and 


figures  2  to  4. 
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REACTION  ON  AN  AIRFOIL  DUE  TO  MOTION  OF  AN  ARBITRARY 
TIME -DEPENDENT  FUNCTION 


In  reference  I  it  has  been  pointed  out  that,  if  the  solution  of  a 
linear,  ordinary  or  partial,  differential  equation  with  a  simple  boundary 
condition  is  known  and  if  the  given  boundary  condition  can  be  obtained 
by  superposition  of  the  simple  boundary  condition,  the  general  solution 
of  the  same  differential  equation  can  then  be  obtained  by  superimposing 
the  elementary  solution  corresponding  to  the  simple  boundary.  In  the 
case  of  the  initial-value  problem,  the  statement  is  also  true.  With  the 
transient  problem,  the  important  relation  of  superposition  is  the  con¬ 
volution  integral.  This  integral  relation  is  associated  with  many  names 
in  modern  mathematical  physics.  In  the  transient  problem  of  heat  transfer, 
it  is  called  Duhamel' s  integral  in  honor  of  the  French  mathematician 
Duhamel  (reference  7,  pp.  ^3-^'^) .  In  England  it  is  commonly  known  as 
Borel's  theorem  (reference  8,  pp.  321-328).  In  this  country  it  is 
commonly  known  as  Carson' s  integral,  particularly  in  electrical  trans¬ 
mission.  In  Germany  it  is  usually  called  the  Faltung  theorem  (refer¬ 
ence  9,  pp.  159-I67),  and  this  was  known  to  Tricomi.  The  concept  of  the 
convolution  Integral  plays  a  very  important  part  in  applied  mathematics. 

In  aerodynamics,  Jones  applied  this  integral  to  the  problem  of  air¬ 
plane  dynamics  in  Incompressible  flow  (reference  10).  The  results  may 
be  considered  natural  in  his  problem  of  ordinary  differential  equations. 
Garrick  (reference  11)  applied  this  integral  to  find  the  relation  of  the 
Wagner  Function  to  the  Theodorsen  function  in  the  case  of  two-dimensional 
incompressible  nonstationary  flow. 

In  reference  1,  the  present  author  constructed  an  integral  relation 
entirely  from  the  physical  concept  because  in  this  problem  of  supersonic 
transient  flow  the  kernel  function  is  discontinuous  in  the  first  deriva¬ 
tive  at  a  certain  Instant  of  time.  The  ordinary  concept  of  the  convolu¬ 
tion  integral  cannot  be  applied  without  detailed  examination,  although 
the  result  can  be  rewritten  in  quite  a  similar  form  to  the  ordinary 
convolution  integral.  For  example,  the  pressure  coefficient  Cp(t’,x'), 

the  lift  coefficient  CL(t' ) ,  and  the  moment  coefficient  C^(t' )  for 

/  \  h(t)  g(t) 

any  arbitrary  a(t)  =  -  or  -  can  be  obtained  by  means  of  the 

U  U 

convolution  integral  of  the  kernel  fxmction  due  to  unit-step  change  in 
angle  of  attack  as  follows: 


4 
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The  subscript  i  should  he  equal  to  1,  2^  or  3  according  to  the 
above  range  of  time  in  the  term  of  the  left-hand  side.  Owing  to  the 
discontinuity  of  these  functions,  the  integral  on  the  right  must  be 
broken  up  into  parts  according  to  the  range  of  time  for  each  fiinction. 
For  example,  in  the  flapping  case,  the  lift  is; 

(a)  For  ^  f, 

M  “  1 


~  %( -0^  Fl(t')  +  I  aj^(T' )F2(t'  -  T' )  dr*  + 

Jo 


(16) 


12 


MCA  TN  2333 


("b)  For  ^  t'  ^ 

M  +  1  M  -  1 


CL(t’ )ii 


^^{+0)  -  ai^(-0y|F2(t'  )  + 


)F2(t' 


)  dr'  + 


a^{T')F^iV  -  t')  dr 


(c)  For  0  ^  t'  ^  M"+T^ 


CtU’) 


=  [^(+0)  -  a^(-0^F3(t’)  +  I 

Jo 


^(T')F3(f 


)  dr' 


Tables  2  and  3  give  the  expressions  for  f ’ s^  F’s,  and  G' s  in  the 
case  of  flapping  and  vertical  gusts.  In  the  case  of  pitching,  the  above 
relation  is  the  same  as  equation  (15),  except  that  a.  should  be  replaced 
by  d  as  given  in  reference  1  except  for  Cp  which  is  replaced  by 

^  in  eqmtion  (15).  Another  complication  in  this  pitching  case  is  the 

x' 

location  of  the  axis  of  rotation  Xq',  which  is  the  ratio  of  the  axis 
location  from  the  leading  edge  divided  by  the  chord.  The  contribution 
due  to  the  nonzero  value  of  Xq'  is  denoted  by  ^F,^,  and  AG^ 

and  should  be  added  respectively  to  f^,  F^,  and  G^.  As  expected, 

=  f,  ,  =  Fv,,  and  =  Gv,  as  shown  in  figures  2,  3)  and  h, 

Xq  Xq  AQ 

respectively. 

As  an  application  of  the  above  convolution  integrals,  the  cases  of 
the  harmonically  oscillating  airfoil  starting  abruptly  from  rest  at  t  =  0 
have  been  investigated  as  follows; 
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Pitching  oscillation  about  the  leading  edge  (xq'  =  O).-  In  the  case 
of  pitching  oscillation  about  the  leading  edge, 

a(  t ' )  =  a 

a(t’)  =  Sel^  = 

CD  -f-TT 

where  —  is  the  frequency  of  oscillation,  3  =  — ,  and  t'  =  — .  (See 
2it  U  C 

table  4.)  Actually,  equation  (19)  is  a  complex  quantity,  and  the  real 
part  is 

Real  a(  t  * )  =  a 
Real  a.(t'  )  =  "a  cos  31' 


The  expression  a(t' )  =  a  for  t'  <  0  is  necessary  in  order  that 
d(t')  be  finite  at  t'  =  0,  because  only  finite  d(t')  is  allowed  for 
pitching.  (For  details  see  reference  1.)  The  imaginary  part  is 


Imaginary  a(t')  =  0  for  t'  <  0 

Imaginary  a(t' )  =  a  sin  3t'  for  t'  >  0 


(21) 


With  the  complex  a(t'),  the  analysis  is  more  convenient  than  with 
the  real  or  the  imaginary  part  alone.  If  the  axis  of  rotation  is  not 
at  the  leading  edge  Xq'  ^0,  the  additional  effect  can  be  obtained  from 
the  flapping  case. 

Flapping  oscillation.-  In  the  case  of  flapping  oscillation. 
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ik 


The  solution  is  given  in  table  5  for  the  three  time  zones.  Some  of  the 
integrals  are  given  in  appendix  B. 

Harmonically  oscillating  gust.-  In  the  case  of  a  harmonically 
oscillating  gust^ 


o-gCt)  =0  t  <  0 

ag(t)  =  t  >  0 


(23) 


The  solution  Is  given  in  table  6.  Some  of  the  integrals  are  given 
in  appendix  B.  All  three  cases  are  essential  to  the  supersonic  flutter. 
The  present  analysis  gives  the  aerodynamic  behavior  of  the  airfoil  for 
the  whole  time  history,  if  the  harmonic  oscillation  in  flapping  and 

M 

pitching  start  abruptly  at  t  =  0.  If  t'  >  —  — the  present  work 
should  check  with  reference  6  exactly.  Appendix  C  shows  the  comparison. 

To  show  the  present  analysis  graphically^  and  have  been 

calculated  with  M  =  1.5  and  p  =  it/2  for  the  three  oscillations. 
Figures  11  and  12  show  the  and  Cj^  of  the  pitching  oscillation. 

Figures  13  and  ik  show  them  for  flapping  and  figures  15  and  l6  show  them 
for  vertical  gust.  The  corresponding  values  of  a(t*)  are  also  shown 
for  comparison  of  phase  shifts.  The  maximum  and  Cjyj  are  larger 

in  the  transient  beginning  than  the  steady  case  which  is  represented  by 

dotted  curves.  The  transient  effect  dies  out  completely  when  t*  >  — ^ — . 

M  —  1 

It  is  interesting  to  note  that  the  extremes  (maximum  or  minimum)  of 
or  always  lead  the  corresponding  extremes  of  a.  But  in  the  case 

of  flapping  oscillation  and  oscillating  gust,  the  extremes  of  or 

always  lag  behind  the  corresponding  extremes  of  a.  As  the  supersonic 
Mach  number  nears  1,  the  transient  effect  becomes  more  pronounced  and 
lasts  longer. 

In  the  time  zone  III  ^0  <  t'  ^  ^  ,  the  expressions  for  Cj^(t*) 

and  )  are  very  simple  in  all  three  cases  of  harmonic  oscillation 

for  pitching,  flapping,  and  vertical  gust.  (Refer  to  tables  4,  5^  and  6.) 

In  the  time  zone  II  ( — ^  t'  ^  — — — j,  the  expressions  are  rather 

Vm  +  1  M  «  1/" 

complicated  and  cannot  be  represented  in  a  closed  form.  Three  new 
functions  or  integrals  have  to  be  defined.  Let 
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m2  -  1 

where  m  =  ^ 6  >  0, 

M  ^  ^ 

function  of  time  t' .  As  designated  ty  Garrick  and  Rubinow, 

may  be  called  the  incomplete  Bessel  function  of  the  zero  order.  It  is 
complex  except  when  =  0.  At  63^  =  0,  eqiiation  (24)  reduces  to  the 

ordinary  Bessel  function  of  the  zero  order  Jo(m)  real.  Simi¬ 

larly^  the  incomplete  Bessel  function  of  the  first  order  can  be  defined 
as 


which  is  also  complex,  but  reduces  to  "the  ordinary  Bessel 

function  of  the  first  order  when  Sp  =  0.  The  new  integral  is  now 
defined  as 


C(P,M;  ep) 


When  @2.  =  it  reduces  to 


(26) 


C(3,M)  =1  de  exp(j^.^i^)  (27) 

For  convenience,  C(p,M)  is  called  the  C-fiinction  and  C(p,M;  dp) 

is  called  the  incomplete  C-function.  The  incomplete  C-function  occurs 
in  time  zone  II  for  all  three  cases.  In  time  zone  III,  it  reduces  to 
C(p,M).  This  function,  to  the  author’s  knowledge,  has  never  been 
explored  before  and  is  investigated  in  appendix  D.  In  time  zone  III, 
Cp_(t’)  and  Cj^(t'  )  can  be  expressed  with  the  known  fiinctions 
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except  C(P^M)  in  all  the  three  cases.  It  is  expected  that  C(p^M) 
is  as  important  to  supersonic  flutter  as  the  Theodorsen  function  is  to 
flutter  in  incompressible  flow. 


TRANSIENT  AERODYNAMIC  BEHAVIOR  OF  THE  CONTROL  SURFACE 


Since  the  principle  of  superposition  holds  for  the  two-dimensional 
linear  problem  if  the  deflection  angle  &  of  the  control  surface  is 
measured  from  the  main  wing,  the  control  surface  itself  may  be  considered 
as  an  independent  airfoil  at  the  corresponding  deflection  angle  6  and 
with  its  chord  C3.  Under  such  a  consideration,  the  lift  and  moment  of 
the  control  surface  itself  in  nonstationary  motion  can  be  obtained 
directly  from  the  result  of  the  airfoil  in  the  early  section,  if  the 
proper  time  scale  is  used.  The  time  required  to  travel  a  length  Cg 

is  t*  =  and,  if  t*  is  used  as  the  time  unit,  the  nondimens ional 

1  connected  with  the  true  time  and  nondimens  ional 

the  relation 


(t"  -  to")t*  =  t  =  t’t 

and 


^11  t  Ut  . 
time  t  =  —  =  —  i£ 

t*  Co 


4.,  t  ut  . 
time  t’  =  rr  =  —  in 

t  C 


(28) 


where  ^9"^*  difference  in  starting  time  of  the  control  movement 

from  that  of  the  main  wing.  Thus  for  an  individual  increment  in  6  as 
the  step  function  at  time  ti”,  the  lift  and  moment  coefficients  of  the 
control  surface  at  later  times  are 


CL(t")  =  A6i(ti")Fj(t"  -  ti") 

=  A5i(ti")Gj(t"  -  ti") 


(j  =  1,2,3) 


(29) 


where  Fj  and  Gj  are  the  same  as  shown  in  table  1  except  that  t*  is 
replaced  by  t". 
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The  lift  coefficient  of  the  control  surface  alone  is 


•t" 


cjt”) 


i£^{o]w  At") 


+ 


(30) 


Then^  the  increment  of  the  total  ving  lift  coefficient  due  to  the  control 
surface  is 
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Similarly^  the  increment  of  total  moment  coefficient  due  to  the 
control  surface  can  he  found* 

As  an  example  of  the  effect  of  the  control  siirface  on  the  total 
lift^  the  following  calculation  is  made. 

(1)  M  =  1.5  a  =  1000  ft/sec  t^  =  0.02  sec  =  5) 

C  =  6  ft  Cs  =  2  ft  (kc  =  1/3)  Xq  =  0  X;j^  =  !(•  ft 

(2)  Operation  schedule  of  the  main  wing: 

a  =  a  =  0  t*<0 

a  =  d  O^t’^t]^* 

d  =  0  t^’  <  t’ 

(3)  Operation  schedule  of  the  control  surface: 


S(t‘) 

=  6(t')  =  0 

t’  <  0 

&(t’) 

=  2S 

ti' 

0  <  t'  ^ 

-  3 

5(t') 

=  0 

U'  <^.  <2tl* 
3  "  3 

6(t') 

=  -2§, 

— <  t'  $  tn  ' 

3  ^ 

6(t' ) 

=  8(t* )  =  0 

t’  > 

The  distribution  of  d(t*)  and  6(tM  against  t’  and  the  geometry 
of  the  airfoil  and  control  surface  are  given  in  figure  17.  The  contri¬ 
bution  of  the  control  surface  to  total  lift^  that  is_,  ACL(t*)^  is.  shown 

in  the  lower  curve;  CL(t’  )  of  the  main  wing  is  shown  in  the  dotted 

curve  which  has  been  given  in  reference  1.  Since  the  ratio  of 

and  d(t*)  is  given^  this  curve  can  be  used  for  any  arbitrary  &.  In 

^tf 

(m2  -  1)^/^ 


this  curve  C 


is  used  instead  of  d. 
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DISCUSSION 


For  the  airfoil  in  a  flow  of  constant  supersonic  speed,  the  tran¬ 
sient  effect  due  to  pitching,  flapping,  or  vertical  gust  will  damp  out 


in  a  time  period  - 

U  -  a 


immediately  after  the  change  in  angle  of  attack 


ceases.  Although  the  result  is  obtained  from  the  linearized  theory,  it 
is  expected  to  be  approximately  true  for  the  nonlinear  theory  if  the 
angle  of  attack  is  reasonably  small.  After  that  time  period,  the  lift, 
wave  drag,  and  moment  become  the  same  as  given  by  Ackeret  in  steady  two- 
dimensional  linearized  supersonic  flow.  The  transient  effect  on  Cp, 

C^,  and  Cm  becomes  more  pronounced  and  lasts  longer  as  the  supersonic 


Mach  number  approaches  1.  The  same  is  also  true  of  the  transient  period 
of  the  harmonic  oscillations  in  pitching^  flapping,  or  vertical  gust 
which  start  from  rest  at  t  =  0. 


In  the  case  of  pitching  with  constant  rate. 


'pa 


'La 


and 


-  dQ^ 


(di  -  do) 


can  be  approximated  satisfactorily  with  a  straight  line 


in  the  time  zone  II,  if  M  >  1.5.  When  such  an  approximation  is  adopted, 
the  case  of  pitching  harmonic  oscillation  can  be  evaluated  very  quickly. 
It  is  also  easy  to  evaluate  any  arbitrary  motion  in  pitching.  In  the 
present  analysis,  the  exact  expressions  in  table  3  were  used  instead  of 
the  above  approximation  because  no  simple  approximation  can  be  obtained 
for  either  flapping  or  vertical  gust. 


In  solving  the  transient  problems  with  the  convolution  integral, 
one  new  function  C(p,M)  is  discovered.  For  the  time  from  the  abrupt 

(M  C\  ' 

f  >  jjj - ^  Cp,  C^,  and  in  the  cases  of  harmonic  oscil¬ 

lations  have  to  be  expressed  in  terms  of  d(p,M).  It  seems  of  comparable 
importance  to  supersonic  flutter  as  the  Theodorsen  function  is  to  flutter 
in  incompressible  flow.  In  the  transient  time  zone  II,  where 

^  t  ^  a  new  function  c(p,M;  ©p) ,  called  the  Incomplete 

M+IU  M-IU  ^  ' 

C-f unction,  occurs  which  assumes  an  importance  equal  to  that  of  the 
incomplete  Bessel  functions.  More  complete  calculation  of  C(3,M)  seems 
useful  for  the  analysis  of ' supersonic  flutter  in  order  to  cover  wider 
ranges  of  Mach  number  and  frequency  than  those  given  in  reference  6. 
Appendix  D  gives  some  of  the  properties  of  C(3,M),  C(p,M;  and 

the  incomplete  Bessel  functions. 
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The  pressiare  coefficients 
Tg 


— rr - ^  in  pitching,  in 

X  (tti  -  ao) 


flapping,  and  ;;; —  in  vertical  gust  are  functions  of  Mach  numher  and 

ttg 

»  t'  -  tn' 

a  conical  parameter  - —  --  only.  They  are  analogous  to  the  behavior 


of  Busemann' s  conical  flow.  As  pointed  out  in  reference  1,  these  results 
can  he  applied  to  a  yawing  infinite  wing,  if  the  leading  edge  is  ahead 
of  the  Mach  line. 


The  effects  of  additional  degrees  of  freedom  such  as  the  control 
surface  or  servoflaps  can  he  evaluated  with  the  result  of  the  pitching 
and  flapping  of  the  main  wing  as  shown  under  Transient  Aerodynamic 
Behavior  of  the  Control  Surface.  With  the  present  basic  approach,  a 
good  aeronautical  engineer  with  ingenuity  should  he  able  to  solve  all 
two-dimensional  problems  of  flutter  and  any  other  arbitrary  motion. 

Of  course  in  the  case  of  complicated  time-dependent  functions  of  angle 
of  attack  some  numerical  or  graphical  Integration  of  the  convolution 
integral  might  be  necessary  as  shown  in  references  1  and  10. 

For  the  case  of  the  harmonic  oscillation  with  constant  maximum 
amplitude  with  abrupt  start  at  t  =  0  the  absolute  magnitudes  of  the 
extremes  in  C^  and  %  for  pitching  and  flapping  degrees  of  freedom 

are  larger  in  the  transient  region  — S  t  S  — —  -  than  those 

M+IU  M-IU 

M  C 

in  the  steady  region  - - -  —  <  t^  particularly  as  the  supersonic  Mach 

M  «  1  U 

number  nears  1.  In  addition^  the  maximum  and  minimum  of  the  oscillating 
angle  of  attack  are  not  in  phase  with  the  maximiim  and  minim\im  of 
or  Cjyj.  In  the  case  of  pitching^  the  load  leads  the  angle  of  attack; 
in  the  case  of  flapping^  the  angle  of  attack  leads  the  load.  Also^  the 
angle  of  attack  leads  the  load  in  the  case  of  vertical  gust. 

As  an  interesting  example  of  the  harmonic  oscillation  building  up 
to  flutter  or  damping  out^  the  case  with  complex  cd  should  be  inves¬ 
tigated.  With  the  convolution  integral,  it  seems  within  the  reach  of 
the  present  analysis.  Of  course,  the  imaginary  part  of  o)  must  be 
determined  from  the  interaction  of  aerodynamic  forces  and  the  elastic 
behavior  of  the  wing. 


The  Johns  Hopkins  University 

Baltimore,  Md. ,  November  4,  19^9 
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3 

6 


T) 

e 


01 


X  = 


frequency  parameter  for  oscillating  airfoils 
deflection  angle  of  control  surface 
source  location  along  y-axis 
running  variable 
lower  limit  of  6 

-1) 


V  x/u(t  -  to)=  x'/(t’  -  to’) 

5  source  location  along  x~axis 

a  =  pM 


T 


0 

00 

Subscripts: 


time  interval 

nondimens ional  time  interval  ('T'/t) 
velocity  potential 

frequency  parameter  (angular  velocity  per  second) 


g  vertical  gust 

h  flapping 

a  pitching  (changing  angle  of  attack) 


6 


control  surface 


2k 


MCA  TN  2333 


APPENDIX  B 

SUMMARY  OF  INTEGRALS 

Important  Integrals  in  the  Case  of  Steady  Harmonic  Oscillation 


The  following  integrals  are  used  in  evaluating  C^  and  Cj^  for  the 

time  zone  t’  > — - —  (as  shown  in  tables  k.  5.  and  6). 

M  ~  1  " 


M-l 


'f  - 


M 

M 

M-1 


’  )  |l  +  “  sin"^  [m(t’  -  T*  )  - 


dr 


-2 


iP 


(i3)‘ 


4 


eiP-"'  }/(t’  -  t'  )2 


}f(t^  -  t'  -  1)^  dj'  = 


Mjtl 

ip)|M2  -  1 


-  M2(t’  -  t’  -  1)^  dT* 


Mjti 


iP 


(*'-  i) 


e  -  f)'!© 


M-1 


eipt* 

ip 


(  PM 
\M  -  cc 


cos  9 


) 


de 


4 
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Important  Integrals  In  the  Case  of  Transient 
Harmonic  Oscillation 

The  following  integrals  are  used  in  calculating  Cj^  and  Cjj  for 

the  time  zone  — ^  t’  ^  —  (as  shown  in  tables  k,  5,  and  6). 

M  +  1  M  -  1 


4- 
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« 


M+l 


ipT 


(f  -  T’ )  \|(t’  -  T’  )^  “  M^(t’  -  t'  -  1)2  dr'  =  - - - - Jl  -  (mf  - 


M)^- 


M  e 


ip 


ip 


nit 


M  cos  0  2  cos  0\  -ipSPg  ^ 


Ucos"-^  •  -M) 


-  cos 

ip 


-1  M(t»  > 
t‘ 


1) 


Pt’  - 


I6' 


e 


ip^'(t'  -  T-)  cos-^  M  1  dr’  =iL_i5_e^'^ 


(‘■-&)  V 


t'  -  t' 


UO 


i3  M  +  1 


iP 


cos 


M(t^  -  1) 


t' 


,  A'-i) 


■  nCOS 

■iP-lfT^  1 


. .  .  I  e  "  “  d0  +  ^  I.' 

-  1  Jcos-l(mf-M) 


I 


t 


7 


M 

M+l 


t' cos“^ 


t'  -  -  1 

t’  -  t' 


dr 


e 


*.(*■-»)  f 


(«2  -  i) 


3/2 


cos"^  ( mt'  -M) 


COS  0 

— i-B  —  p 

(M  +  cos  0)e  ^  Tp^  ^6* 


28 


MCA  TN  2333 


* 


4 


+  1  +  2M  cos  <l>  + 


^  cos  A 
13  7 


c^ 

!  ^  ^  m  dj!!  + 


( mt' -M) 


MCA  TN  2333 


29 
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APPENDIX  C 


► 

CHECK  WITH  MCA  TN  II58 


To  check  the  present  results  through  convolution  integrals  with 
Garrick  and  Ruhinow's  work,  a  case  of  the  steady  harmonic  oscillation 
Of  the  angle  of  attack  about  the  leading  edge  is  investigated.  With 
their  notation,  Cj^  due  to  the  harmonic  oscillation  of  the  angle  of 

attack  alone  (a  =  a^e^^j  may  be  verified  from  their  equation  (26)  as 


^La  “  (Cl) 

At  Xq  =  0, 

L3  =  L3'  L3  =  13'  -  2xoLi 

L4  =  =  L]^'  -  2xqL2 


At  M  -  2,  for  i  =  1,667,  from  their  table  II, 


L3'  =  1.50219  Lj^'  =  0.69968 


Hence, 


Cla  =  0^(2.16315  +  1.00T54i)ei^ 
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NOTATION 


whence,  at  M  =  2  and  p  -  1.2,  C^ 
The  expressions  are  identical,  within 


a  exp  ( ipt ’ ) (2. 16315  +  1.007531) 

computational  limits* 
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APPENDIX  D 

INVESTIGATION  OF  THE  FEW  IMPORTANT  INTEGRALS  RELATED  TO 
SUPERSONIC  FLUTTER  AND  TRANSIENT  PROBLEMS 
C -Function 


The  problem  is  investigation  of  the  C-function  or  the  integral 


C(p,M) 


(3  >  0) 
(M  >  1) 


If  the  above  integral  is  differentiated  with  respect  to  3  and  M, 


^3^ 


a  r  dJexp  (  ] 

« Jo  L 


cos  6 


(M  -  cos  6) 


M 


cos  6 


7^  (m"  «)] 


(m  -  cos  e) 


(Dl) 


(D2) 


(D3) 


With  the  above  three  relations,  it  is  found  that  C(3^M)  =  Cp  +  iCj 
satisfies  the  differential  equation 


M  =  1  ^  +  M  ^ 
S3  '  a3^  p  ^ 


(d4) 
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If  the  real  and  imaginary  parts  of  C(3,M)  are  taken  separately, 

M 

3  ^ 

^  +  M  ^ 

Sp  ^p2  p  ^ 

which  are  two  simultaneous  differential  equations.  The  required  six 
boundary  conditions  can  be  obtained  from  C(p,M)  as  follows. 

For  P  =  0, 


For 


Cr(o,m)  =  1 

C];(0,M)  =  0 

il 

HI 

^  1 
ro 

Cp  =  cos  p 

Cj  =  -sin  p 

(d6) 


(D7) 


The  integral 


can  be  written  as 


where  a  =  PM.  Under  the  new  definition,  the  function  C(or,M)  can  be 
shown  to  satisfy  the  differential  equation 


(D8)  f 
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SCj 


It  is  easy  to  show  that  the  real  and  imaginary  parts  of  C(a^M) 
must  satisfy  the  equations 


^  L 


a  ^  a  0  ^  a  ^ 

Scr^  \ 


^2  /  a2Cj\ 

a^a-^  =_^  l^M<oo 

Bff  \  ha^.  /  hit 


The  corresponding  boundary  conditions  are: 
For  CT  =  0,  M  >  1, 


Cr  =  1 


Cj  =  0 


|m^  -  1 


(Dll) 


(m2  -  l' 


For  M  ^  00^ 


Cr  =  1 


Cj  =  0 


(D12) 
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3^ 

The  preceding  function  is  so  important  to  the  supersonic  flutter 
problem  that  thorough  investigation  of  the  function  and  the  other  ^ 

solutions  of  the  corresponding  differential  equation  seems  advisable. 

To  show  the  general  behavior  of  this  new  function^  a  case  with  M  =  1.3 
and  0  <  ^  %  3  has  been  calculated  numerically.  See  table  7  and 
figure  l8. 

Incomplete  C-function 
The  incomplete  C-fimction  is 


C(P,M;  ei) 


(D13) 


It  can  be  shown  that  this  incomplete  C-function  also  satisfies  the 
differential  equations  (l)5)  except  the  boundary  conditions  are  obtained 
from  equation  (D13). 

For  3=0^ 


(D15) 


To  show  the  nature  of  the  incomplete  C-function,  a  case  with  M  =  1.5, 
p  =  and  0  jt  has  heen  calculated  hy  means  of  numerical  inte¬ 

gration.  See  table  8  and  figure  19.  The  above  data  are  used  in  the  tran¬ 
sient  behavior  as  shown  in  figures  11  to  l6. 
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Incomplete  Bessel  F-unctions 
The  incomplete  Bessel  functions  are 


Ojt 

^  cos  0)  d0  (DI6) 

J^i 


pjr 

Jl(z^0l)  =  cos  0  exp  (+iz  cos  0)  (10  (DI7) 

J®1 


These  two  functions  become  Jq(z)  and  J^(z),  respectively,  when 
0]^  =  0.  They  have  been  called  the  incomplete  Bessel  functions  by 
Garrick  and  Rubinow,  The  case  M  =  I.5,  P  =  f ^  0  ^  03_  ^  jt  has 

been  calculated  by  numerical  methods.  It  is  interesting  that  each  has 
a  real  and  an  imaginary  part.  At  0^  =  0,  both  become  real  and  ordinary 
Bessel  functions.  Both  play  as  important  roles  as  the  incomplete 
C-function  in  harmonic  oscillations  at  the  transient  time  zone  II.  See 
table  8  and  figures  20  and’  21. 


36 


NACA  TN  2333 


REFERENCES 


1.  Chang,  Chieh-Chien:  The  Transient  Reaction  of  an  Airfoil  Due  to 

Change  in  Angle  of  Attack  at  Supersonic  Speed.  Jour.  Aero.  Sci., 
vol.  15,  no.  11,  Nov.  1948,  pp.  635-655. 

2.  Biot,  M.  A.:  Loads  on  a  Supersonic  Wing  Striking  a  Sharp-Edged 

Gust.  Rep.  No.  SA-247-S-7,  Cornell  Aero.  Lah,,  Cornell  Res. 

Foundation,  Inc.,  1948. 

3.  Schwarz,  L. :  Plane  Nonstationary  Theory  of  the  Wing  at  Supersonic 

Speed;  an  Abstract  of  a  Report  with  the  Same  Title.  Translation 
No.  F-TS-934-RE  (ATI  22564),  Air  Materiel  Command,  Army  Air  Forces, 

March  1947. 

4.  Heaslet,  Max.  A.,  and  Lomax,  Harvard;  Two-Dimensional  Unsteady  Lift 

Problems  in  Supersonic  Flight.  NACA  Rep.  945,  1949.  (Formerly 
NACA  TN  1621.) 

5.  Miles,  John  W. :  Transient  Loading  of  Airfoils  at  Supersonic  Speeds. 

Jour.  Aero.  Sci.,  vol.  15,  no.  10,  Oct.  1948,  pp.  592-598. 

6.  Garrick,  I.  E.,  and  Rubinow,  S.  I.:  Flutter  and  Oscillating  Air- 

Force  Calculations  for  an  Airfoil  in  a  Two-Dimensional  Supersonic 
Flow.  NACA  TN  1158,  1946. 

7.  V.  Kannan,  Theodore,  and  Blot,  Mairrice  A.:  Mathematical  Methods  in  ,♦ 

Engineering.  McGraw-Hill  Book  Co.,  Inc.,  194o. 

8.  McLachlan,  N.  W. :  Complex  Variable  and  Operational  Calculus  with  * 

Technical  Applications.  The  MacMillan  Co. ,  1946. 

9.  Doetsch,  G. :  Laplace  Transformation.  Dover  Publications,  1943. 

10.  Jones,  Robert  T.:  Calculation  of  the  Motion  of  an  Airplane  -under 
the  Influence  of  Irregular  Disturbances.  Jour.  Aero.  Sci., 
vol.  3,  no.  12,  Oct.  1936,  pp.  419-425;  A  Simplified  Application 
of  the  Method  of  Operators  to  the  Calculation  of  Distiorbed  Motions 
of  an  Airplane.  NACA  Rep.  56O,  1936. 

Garrick,  I.  E.:  On  Some  Reciprocal  Relations  in  the  Theory  of  Non¬ 
stationary  Flows.  NACA  Rep.  629^  1938. 


11. 


MCA  TN  2333 


37 
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CHAEACTERISTICS  OF  A  FLAPPING  AIRFOIL  WITH  K  (t  >  0}  OR  ADDITIONAL 
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TABLE  4 

CHANGING  ANGU:  OF  ATTACK 


El  =  oe^P^’  t'  >  0 

in  the  case  Xq’  =  0.  See  figs.  H  and  12. 

[,  =  a  0  >  t' 


FLAPPING  WING 
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'c-fimctlon  (3,M). 

Incomplete  C-function  C(3^M;  ©i),;  see  appendix  D. 
Incomplete  Bessel  functions;  see  appendix  D, 


VERTICAL  GUST 


k2 
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C -function. 

Tncomplate  C-function  C(3,M;  ©i);  see  appendix  D. 
Incomplete  Bessel  functions;  see  appendix  D, 
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TABLE  7 

VALUES  OF  C-FUNCTION  C(p,M)  FOR  M  =  1.3 


p 

Real 

Imaginary 

0 

1.0 

0 

.3 

.84153 

-.417961 

.5 

. 62101 

-.567971 

.8 

. 32000 

-.569621 

1.0 

.21844 

-.498441 

1-2 

.19147 

-.448541 

1.5 

.17883 

-.477121 

1.7 

.12058 

-.534231 

o 

OJ 

-.48081 

-.562791 

2.2 

-.16005 

-.511541 

2.4 

-.22697 

-.424911 

2.5 

-.23982 

-.381231 

2.93 

-.22265 

-.293551 

kk 
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TABLE  8 

VALUES  OF  INCOMPLETE  C-FUNCTION  AND  INCOMPLETE  BESSEL  FUNCTIONS 

[p  =  n/2,  M  =  1.5] 


Sl  , 

C(P,M 

>  ®l) 

Jo  (5 

>1) 

Jl(mJ®l) 

0 

-0.05339 

-0. 471481 

0.29055 

0 

0,58i48 

0 

.2 

-.04952 

-.534931 

. 30946 

,060781 

.52110 

. 018801 

-.02526 

-.593171 

.32378 

.122781 

.  46198 

.032491 

.  6 

. 02500 

-.630791 

. 32891 

,.  I86l4l 

,40647 

.037071 

.8 

.08770 

-.634931 

. 32058 

,249101 

,35837 

.030881 

1.0 

•  .14434 

-.607201 

.29582 

.307501 

. 32204 

.015731 

1.2 

.18435 

-.558191 

.25417 

.355361 

. 30020 

-.002891 

1.4 

.20575 

-.498481 

.19871 

. 385801 

.29172 

-.015721 

1.6 

.21084 

-.435171 

.13598 

. 394201 

.29074 

-.021901 

1.8 

.20310 

-.372061 

.07454 

.379001 

.28841 

-.014091 

2.0 

. 18569 

-.310871 

.02257 

. 342811 

.27643 

,002481 

2.2 

.16124 

-.252121 

-.01429 

.291251 

,25027 

.020821 

2.4 

.13173 

-.195731 

-.03407 

,230961 

.21010 

.033771 

2.6 

.09867 

-.141331 

-,03808 

.167561 

.15938 

.036831 

2.8 

.06326 

-. o8844i 

-,02980 

, 104501 

.10248 

.029291 

3.0 

.02644 

-.036501 

-.01367 

.042951 

.04279 

.013591 

3.1 

,00778 

-,010711 

-.004o8 

.012591 

,01259 

.004111 

rt 

0 

0 

0 

0 

0 

_ 

0 

_ 
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Figure  4.-  Plot  of  Gu  or  — 2_ 

Xq' 

numbers.  Gj,  =  ; 

“h 


as  a  function  of  t '  -  tQ  '  at  various  Mach 


aG^',  contribution  to 


^Md 

ttf  -  Oq 
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0  2  4  6  8 


t' 

Figure  17."  Distribution  of  d  (f)  and  6(t')  against  t'  and  geometry  of 
airfoil  and  control  surface.  Contribi^ion  of  control  surface  to  total  lift 
shown  by  curve  at  bottom  of  figure.  CL(t')  shown  by  dotted  curve 
(from  reference  1).  t^'  =  5. 
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Figure  19.-  Incomplete  C-function  C(p,M;  9^)  for  p  =  n/2  and  M  =  1.5. 


Figure  21.-  Incomplete  Bessel  function  of  first  order 

M  =  1.5  and  p  =  n/2. 
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